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Abstract 

We suggest a new strategy for proving large N duality by interpreting Gromov- 
Witten, Donaldson-Thomas and Chern-Simons invariants of a Calabi-Yau threefold 
as different characterizations of the same holomorphic function. For the resolved 
conifold this function turns out to be the quantum Barnes function, a natural q- 
deformation of the classical one that in its turn generalizes Euler's gamma function. 
Our reasoning is based on a new formula for this function that expresses it as a 
graded product of g-shifted multifactorials. 

Introduction 

What is the topological string partition function of the resolved conifold? We should 
explain that heuristically one can assign string theories to each Calabi-Yau threefold and 
some of them such as topological A-models [62], only depend on its Kahler structure. Their 
topologically invariant amplitudes are then collected into a generating function called the 
partition function. Remarkably, this partition function may remain unchanged even if a 
threefold undergoes a topology changing transition [2T] . 

A traditional approach is to interpret the string partition function as the Gromov- 
Witten partition function. For the resolved conifold X := 0(—l) © O(-l) it was origi- 
nally computed by Faber-Pandharipande [Ej, see also [61] 

oo 

Z' x (a;q) = l[(l-aq n ) n . (1) 

71=1 

Here a = e _t , q = e lx and t, x are known as the Kahler parameter and the string coupling 
constant respectively. In mathematical terms, they are just formal variables and 

oo oo 



g=0 d=l 



where {l) 9t d is the Gromov-Witten invariant of genus g degree d holomorphic curves in 
the resolved conifold. 

The incompleteness of this answer does not reveal itself until one considers dualities 
that relate Gromov-Witten invariants to other invariants of Calabi-Yau threefolds. One 
may notice that (jSJ) is missing degree zero terms (hence the '). This is not a slip, they can 
not be packaged into a form as nice as ([I]). This was not considered much of a problem 
until Donaldson-Thomas theory f25l WI\ came about since degree zero (constant) maps 
are trivial anyway. But apparently dualities have little tolerance for convenient omissions. 
For Gromov-Witten/Donaldson-Thomas duality to hold ([1]) has to be augmented as 

Z X = Z° X Z' X ^MZ' X , (3) 

where 

oo 

M(q) := -?")-" (4) 

n=l 

is the MacMahon function, classically known as the generating function of plane partitions 
|38j . In all honesty, this is not quite true as \nM(e lx ) has some spurious terms in its 
expansion at x = and only accounts for genus g > 2 terms correctly (see Section [2} . 
Also in Donaldson-Thomas theory one has Z X T = M 2 Z' X DT , not Z x = MZ X . In a 
recent reformulation of the Donaldson-Thomas theory the reduced partition function 
Z' DT is even defined directly and the Mac-Mahon function is banished altogether. Let 
us disregard this minor discrepancy for now since even answer ([3]) is incomplete. 

This becomes apparent in light of another duality of Calabi-Yau threefolds, large N 
duality. This one relates Gromov-Witten invariants of the resolved conifold to Chern- 
Simons invariants of the 3-sphere. The usual formulation defines Chern-Simons theory 
as a gauge theory on a Un or SUn bundle over a real 3-manifold M. Less recognized 
despite Witten's famous paper [62] is the fact that it also gives invariants of Calabi- 
Yau threefolds. As Witten pointed out it can be viewed as a theory of open strings 
(holomorphic instantons at oo in his terminology) in the cotangent bundle T*M ending 
on its zero section. T*M is canonically a symplectic manifold (even Kahler if M is 
real- analytic) with first Chern class C\(T*M) = 0, i.e. a Calabi-Yau. In particular, 
T*S 3 is diffeomorphic to a quadric x 2 + y 2 + z 2 + w 2 = 1 in C 4 . One of the reasons 
this interpretation did not get much currency is that the strings in question are very 
degenerate, they are represented by ribbon graphs, not honest holomorphic curves. In 
fact, there are no honest holomorphic curves in T*M at all except for the constant 
ones [331 Another reason perhaps, is that open Gromov-Witten theory is still in 
its infancy and the powerful algebro-geometric techniques that dominate the field can 
not be directly applied. There are successful approaches that replace open invariants 
with relative ones [TH [36] but only as a tool for computing closed invariants. In the 
other direction, there exists a detailed if only formal correspondence between geometry 
of real oriented 3-manifolds and Calabi-Yau threefolds and Donaldson-Thomas theory 
can be seen as a 'holomorphization' of Chern-Simons theory under this correspondence 
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J. Thus, comparing the Chern-Simons partition function Z s -a to Zx promises some 
useful insight. 

Once again ignoring some irrelevant prefactors, Z$3 can be written as Z$z « £~ z Zx, 
where z = itx~ x so that a = q z and 

£(q):=H(l-q n r 1 (5) 

n=l 

is the classical Euler generating function of ordinary partitions. At this point it is ap- 
propriate to introduce notation that allows one to write Z' x , M and £ uniformly. Let 

00 

(o;g)2>:=l- a> (a; q) g := J] (1 - aq^~ +i *) (6) 

be the q-multifactorials then (see Section [5]) 

Z' x (a;q) = (aq;q) ( £, M(q) = £{q) 



00 



Using g and z as variables we see that 

1 



(2) 



(7) 



{q; q)% 

After some thought one may sense a pattern here. We shall see in Section [5] that it makes 
sense to join one more factor to the product and consider 

G q (z + 1) := -^ir (q; q)S z (<f +1 ; g)g. (8) 



This is the quantum Barnes function of Nishizawa jl5] and our candidate for the 
partition function of the resolved conifold. All factors above are required to make it 
transform as 

G q (z + l) = T q (z)G q (z), (9) 

where T q is Jackson's quantum gamma function deforming the classical one. This in 
turn satifies T q (z + 1) = (z) q T q (z) with the so-called quantum number (z) q := -j^-- 
This makes G q a deformation of the classical Barnes function that satifies (Q with q-s 
removed. 

The picture above is cute but not quite true and clear-cut identities are spoiled 
by pesky disturbances discussed in Sections |2] and HI These disturbances are a large part 
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of the reason why large iV duality is so hard to prove even in simple cases. Still, G q 
emerges as a common factor in Gromov-Witten, Donaldson-Thomas and Chern-Simons 
theories (Theorem [3]). One may notice that we conspicuously omitted the most famous 
of Calabi-Yau dualities, mirror symmetry. This is partly because local mirror symmetry 
is poorly developed, and partly because to the extent that its predictions can be divined 
[16] they match the Gromov-Witten ones completely. There is a structural prediction of 
mirror symmetry that seems relevant. For compact Calabi-Yau threefolds Z is predicted 
to have modular properties [12], i.e. obey transformation laws under z i— > z + 1 and 
z i — > —1/z. For open threefolds like the resolved conifold only the first one survives and 
is expressed by (EJ). 

What are we to make of the above chain of augmentations? Perhaps, string theories 
on Calabi-Yau threefolds are only partial reflections of some hidden master-theory. Wit- 
ten's candidate for such a theory is the mysterious M-theory living on a seven- dimensional 
manifold with G2 holonomy that projects to various string theories on Calabi-Yau three- 
folds. Another unifying view of Gromov-Witten and Donaldson-Thomas theories, via 
non-commutative geometry, also emerged recently (56]. Different projections are equiva- 
lent even though they may live on topologically distinct threefolds and reflect the original 
each in their own way. So far we ignored these ways relying instead on magical changes 
of variables. It is time to dwell upon them a bit. This will also serve as our justification 
for spending so much ink on the resolved conifold. 

The relation between Gromov-Witten and Donaldson-Thomas invariants is very sim- 
ple [12]. For the resolved conifold we have 

00 00 

n=0 d=l 

with Z' x the same as in ([2]) and D n d the Donaldson-Thomas invariants. In other words, in 
each degree (— l) n D nt d are simply the Taylor coefficients of Z' x at q = while {l) g ,d are the 
Laurent coefficients in x corresponding to q = 1 with q = e lx . The relation with Chern- 
Simons invariants is more complicated. Traditionally, one has to take q = e 2m /( k + N )^ 
where k, N are the two parameters of Chern-Simons theory, rank and level. They are 
positive integers making q a root of unity. Not all roots of unity are covered in this way 
but more sophisticated formulations allow one to include any root of unity. Naively, if 
the duality conjectures hold Donaldson-Thomas invariants give us an expansion at q = 0, 
Gromov-Witten invariants at q — 1 and Chern-Simons invariants give values at roots of 
unity of more or less the same function. But only naively. First of all, Donaldson-Thomas 
generating functions are a priori only formal power series and may not have a positive 
radius of convergence. We need it to be at least 1 to make a comparison. Things are 
nice in higher degree [51] but in degree zero it is exactly 1 and every point of the unit 
circle is a singularity. This is remedied easily enough in the Gromov-Witten context since 
we can interpret Z x (e lx ) = Yl^Lo^X) gfl x29 ~ 2 as an asymptotic expansion at the natural 
boundary (Section . But Chern-Simons invariants are not graded by degree and the 
degree zero speck turns into a wooden beam spoiling the whole partition function that we 
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wish to evaluate. With the resolved conifold being the simplest non-trivial case we get a 
preview of the difficulties that will arise in general. This brings us to a paradox: for large 
N duality to even make sense the formal power series better converge to holomorphic 
functions extending to the unit circle or at least to roots of unity. This is not the case 
already for Zx and an additional factor in Zg3 appearing in (j7]) is needed to make it 
happen (see comments after Corollary [2]) . This is another reason to accept the quantum 
Barnes function as the completed partition function. 

Since conjecturally Z Y = Ms^' for any Calabi-Yau threefold Y [12] this phenomenon 
is likely to be general. The above discussion suggests that the master-invariant that 
manifests itself through dualities is a holomorphic function on the unit disk. The three 
theories we discussed showcase three different ways to package information about it. The 
dualities reduce to repackaging prescriptions. Physicists developed resummation tech- 
niques that transform generating functions one into another but they lead to unwieldy 
computations for the resolved conifold and do not produce conclusive results even for 
its cyclic quotients [2J. Since repackaging involves transcendental substitutions, analytic 
continuation and asymptotic expansions - things one does with functions and not with 
formal series - it makes sense to identify the underlying holomorphic functions to estab- 
lish a duality. This is the strategy of this paper and it distinguishes it from previous 
approaches [H [211 HI] that use double expansions in genus and degree. This makes for 
a cleaner comparison of partition functions with a clear view of what matches and what 
does not match in them (Theorem [3]). It is also hoped that the idea generalizes to other 
threefolds. 

The paper is organized as follows. Section [T] is a review of basic notions of Gromov- 
Witten theory with emphasis on generating functions. In particular, we note that free 
energy is a shorthand for the Gromov-Witten potential restricted to divisor invariants. 
The well-known irregularities in degree zero are then naturally explained. In Section [2]the 
MacMahon function is examined in detail to determine to what extent it can be viewed 
as the degree zero partition function of Gromov-Witten invariants. We describe resum- 
mation techniques used by physicists, and then recall an old but little known asymptotic 
for it due to Ramanujan and Wright adapting it to our context. Sections [3] and H] give a 
description of the topological vertex and the Reshetikhin-Turaev calculus, diagrammatic 
models that compute Gromov-Witten and Chern-Simons partition functions respectively. 
Similarities between the two are specifically stressed. Section H] ends by expressing both 
partition functions via the quantum Barnes function (Theorem [3]). Since this function 
and its higher analogs are relatively recent (1995) we give a self-contained exposition of 
their theory in Section [5] different from the author's [15] . In particular, we prove the 
alternating formula (jSJ) that connects G q to the Calabi-Yau partition functions and ap- 
pears to be new (Theorem H]). In Conclusions we point out relations between Calabi-Yau 
dualities and holography and share some thoughts and conjectures inspired by the re- 
solved conifold example. The Appendix lists basic properties of the Stirling polynomials 
needed in Section |5j 

Acknowledgements. The author wishes to thank D.Auckly and D.Karp for sharing 
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their thoughts and notes on Gromov-Witten and Chern-Simons theories, and suggesting 
numerous corrections and improvements to the original draft. I am also indebted to 
M.Marino for his extensive email comments on the physical background of Calabi-Yau 
dualities. His influence can be felt in the overall mindset of this paper. 

1 Generating functions of 
Gromov-Witten invariants 

There are a variety of generating functions appearing in the literature: Gromov-Witten 
potential, prepotential, truncated potential, partition function, free energy, etc. In this 
section we briefly review basic definitions from Gromov-Witten theory and relationships 
among some of the above generating functions. Perhaps, the only unconventional notion 
is that of divisor potential which leads most naturally to the free energy and the partition 
function. 

Stable maps. Let X be a Kahler manifold of complex dimension N. We wish to 
consider holomorphic maps / : £ — > X of Riemann surfaces with n marked points into 
X that realize certain homology class a G H 2 (X, Z). The space of such maps is denoted 
Mg >n (X, a). There is a natural (Gromov) topology on this moduli space but it is not 
compact in it. To get Gromov-Witten invariants we need to integrate over the moduli so 
we have to compactify. The appropriate compactification was discovered by Kontsevich 
and its elements are called stable maps. They are holomorphic maps from prestable 
curves, i.e. connected reduced projective curves with at worst ordinary double points 
(nodes) as singularities. A map is stable if its group of automorphisms is finite, i.e. 
there are only finitely many biholomorphisms a : S — > S satisfying / o a — / and 
a(pi) = pi, where p\ . . . ,p n are the marked points. Intuitively, we allow Riemann surfaces 
to degenerate by collapsing loops into points. Since only genus and 1 curves have 
infinitely many automorphisms (Mobius transformations and translations respectively) 
the stability condition is non-vacuous only for them and only if the map / is trivial, i.e. 
maps everything into a point. It requires then that each genus 0(1) component have 
at least 3(1) special points, nodes or marked points. Under favorable circumstances the 
space of stable maps M g>n (X, a) up to reparametrization is itself a closed Kahler orbifold 
of dimension 

dim£ r M g , n (X, a) = ( Cl (X), a) - (N - 3)(g - 1) + n. (11) 

For instance, this is the case if X = CF N and = 0. Above c x {X) is the first Chern 
class of the tangent bundle and (,) the co homology /homology pairing. The notation 
anticipates that in general the moduli is neither smooth nor has the expected dimension 
so pip is called the virtual dimension. A deep result in Gromov-Witten theory asserts 
that despite the complications there is a cycle of expected dimension [Mg >n (X, a)] vir called 
the virtual fundamental class that one can integrate over. 
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Primary invariants. Presence of marked points allows one to define evaluation maps 

ev, : M g>n (X,a) — X 
f — > f{Pi) 

and pullback cohomology classes 7* from X to M 9tn (X, a). These pullbacks are called 
the primary classes on M g>n (X,a) [I9| 132] . The primary Gromov-Witten invariants are 

(7i---7n),, a := /_ evI(7i)U---Uev;( Tn ), (12) 

where U is the usual cup product and the integral denotes pairing with [Mg >n (X, a)] vir . 
Again under favorable circunstances the primary invariants have an enumerative inter- 
pretation. Namely, (71 • • ■ r y n )g,a is the number of genus g holomorphic curves in a class 
a G H 2 (X, Z) passing through generic representatives of cycles Poincare dual to 71, . . . , j n 
[HI ED]- In general, the enumerative interpretation fails and (71 • • • 7„) 9 . a are only rational 
numbers, this is always the case for Calabi-Yau manifolds. Most of the primary invari- 
ants are zero for dimensional reasons. Indeed, the complex degree of the integrand in 
( fl2l) is |(deg7i + • ■ ■ + deg7„) and for the integral to be non-zero it should be equal to 
the virtual dimension (TTTT) . There are other natural classes on M 9jTl (X, a) that lead to 
more general Gromov-Witten invariants, gravitational descendants and Hodge integrals 
[HI [El [19], but we need not concern ourselves with them here. 

It is convenient to arrange the primary invariants into a generating function [30]. To 
this end we note that they are linear in insertions 7j and we can recover all of them from 
(1) 5)Q and those with insertions chosen from an integral basis hi, . . . , h m in H + (X, Z) : = 
Q) n>0 H n (X, Z). One may worry about torsion but torsion classes are not represented by 
holomorphic curves and can be ignored. Thus any (71 • • • 7n)g,a is a linear combination of 
(h 1 ^ ■ ■ ■ h^ , )g t0l , where the 'powers' pi stand for repeating hi that many times. Introduce 
formal variables ti, . . . ,t m for each element of the basis. Heuristically, they represent 
(minus) Kahler volumes of h\,...,h m and are called Kahler parameters, especially in 
physics literature. Analogously, let £1, . . . , ^ be a linear basis in H 2 (X, Z) and Q 1 , . . . , Q k 
the corresponding formal variables. We write {h^ 1 ■ ■ ■ h^) j* with d := (di, . . . , c4) for 
short when a = di£i + • • • + c46s- The numbers di, . . . dk are called degrees. Finally, 
we need one more variable x, the string coupling constant, to incorporate genus. The 
primary Gromov- Witten potential (relative to the above bases choices) is 

~ ~ t Pl ■■■ fP™ 
F(ti,...,t m ;Qi,...,Q k ;x):=J2 E ( h i ' ' ' h ™ >*f in ^ ' ■ ■ Qt ^ ■ 

g=0 Pi,..., P m=o Pl-"-Pm- 

d lt -,d k =0 

(13) 

This particular choice of a generating function is by no means obvious and is inspired 
by two-dimensional topological quantum field theory. The power 2g — 2 instead of just 
g has in mind the Euler characteristic —(2g — 2) of a genus g Riemann surface. For X 
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Kahler T is defined as at least a formal power series in Q[[£j, Qi, x]] [2U]. Under a change 
of bases (h^ 1 ■ ■ ■ h^) g j transforms as a tensor. One may entertain oneself by writing a 
tensor potential that is an invariant, see [30]. In [191 [32] a more general Gromov-Witten 
potential is considered that incorporates gravitational descendants and accordingly has 
more formal variables. 

Let X := 0{— 1) ffi @{— 1) be the resolved conifold, the sum of two tautological 
line bundles over CP 1 . Being a vector bundle over CP 1 it is homotopic to its base 
and has the same homology and cohomology. In particular, H2(X,Z) = Z[CP ] and 
H*(X, Z) = Z[/i]/ (h 2 ), where h is the Poincare dual to the class of a point in CP 1 . Thus, 
H + (X, Z) is spanned by h and H 2 (X, Z) is spanned by £ := [CP 1 ], the fundamental class 
of CP 1 . Hence we need only one t and one Q variable. The primary potential simplifies 
to 

oo 

F(t;Q;x):= (W) g4 - Q d x 2 ^ 2 . (14) 

P ,d,g=o P ' 

Divisor equation and free energy. We will be interested not even in all primary 
invariants but only in those corresponding to combinations of divisor classes, elements 
of H 2 (X,Z). Divisor invariants turn out to be most relevant to large N duality. In 
non-compact manifolds the name is misleading since there is no Poincare duality. For 
example, the hyperplane class of CP 1 is a divisor class in 0(—l) ®0{— 1) despite the fact 
that it is not dual to any divisor. But in closed manifolds divisor classes are precisely 
Poincare duals to divisors, cycles of complex codimension one. Invariants (h^ 1 ■ ■ ■ h^) g g 
containing only divisor classes can be reduced to (l) g g using the so-called divisor 
equation. The latter is one of the universal relations among Gromov-Witten invariants 
coming from universal relations among moduli spaces of stable maps with the same target 
X. One of them is [El [32] 

n*\M g , n -i(X, a)] vh = [M g>n (X, «)] vir , (15) 

where M g , n (X,a) — > M 9) „_i(X, a) is the map forgetting the last marked point. Its 

7T 

consequence is the divisor equation 

(hji ■ ■■Jn)g,a = K a ){ll ' ■•Jn)g,a, (16) 

where h G H 2 (X, Z) and 7« are arbitrary. There are two exceptions to the validity 
of (TT5"j) and hence (TToT) . both in degree zero. If a = then M gtn (X, 0) consists of 
constant maps. The stability condition requires domains of stable maps in this case to 
be themselves stable, not just prestable. But when g = 0(1) a stable curve must have 
at least 3(1) marked points so the spaces of curves M 0)0 , M 0i i, Mo,2, M 1)0 are empty. 
However, M 3 , M^x are n °t and ([15]) fails for (g, n) = (0, 3), (1, 1). 

Since H 2 (X, Z) ~ H2(X, Z) modulo torsion and £i, • • • , £fc form a basis in H2(X, Z) 
there are precisely k basis elements in H 2 (X,1t). We assume without loss of generality 
that hi,...,hk are the ones and that they are dual to £i, . . . , i.e. = <%. The 
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divisor equation may now be used to flush all the insertions out of the divisor invariants. 
By induction from (1161) 

(K---K k ) g j=df...df(l) g ^ (17) 

— * 

assuming d 7^ to avoid low-genus problems in degree zero. Define the truncated divisor 
potential J- d iv (ti, . . . , t&; Q±, . . . , Q k ] x) as in (jl3l) but restricting the sum to pi, . . . ,pk and 
d ^ 0. Using ffTTl) we compute 

Fdivi^ ■ ■ ■ ,tk',Qi, 

00 

9=° d^O 

Obviously, as far as divisor invariants go Qi, . . . , Q k are redundant and we can set them 
equal to 1. This naturally leads to another generating function [26, 42J: 

Definition 1. The reduced Gromov-Witten free energy is 

00 

F'( tl , ...,t k ;x):=J2 E< 1 >^ edlh • • • ^ x29 ~ 2 - ( 18 ) 

9=° d^O 

Its exponent Z'(ti, . . . , t k ; x) := exp(F'(ti, . . . , t k ; x)) is called the reduced Gromov-Witten 
partition function. We write F' x , Z' x when the target manifold needs to be indicated. 

The reduced free energy is non-zero only if (ci(X), a) — (N — 3)(g — 1) = for some 
class a 7^ 0, see (fTTj) . If X is a Calabi-Yau then C\(X) = and if in addition it is 
a threefold then also N = 3 and the non-triviality condition holds for all classes and 
genera. For a toric Calabi-Yau X the reduced partition function Z' x is the quantity 
directly computed by the topological vertex algorithm [HI EU ES, HI] . 

Degree zero. The moduli spaces M g n (X, 0) consist of stable maps mapping stable 
curves into points. Therefore they split [H] 

I s , n (I,0)=I 9 , n xI. (19) 

This reduces degree zero invariants to integrals over the spaces of curves and over X. The 
divisor equation ( fl6l) still holds for n > 4(2) for genus g = 0(1) and for all n in higher 
genus. Moreover, since a = now it directly implies that all the divisor invariants vanish 
except possibly for those that can no longer be reduced. Therefore, in genus g > 2 the only 
surviving invariants are (l) ffl o and in genus 0, 1 we are left with (/j|) ,o, (h"ihj)o,o, (hihjhi)o,o 
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*) = £ £<i>^ ot- Qt E - ltir ' ' ' {dmtmr 



9=0 d^0 



pi,...,p fe =0 



Pl ] - ■•■Pr, 



, Q d k x 2 9 -2 e d ltl . . . e d ktk = ^ J2( 1 ) 9 , d ;(Qi etl ) dl ■ ■ ■ iQke tk ) dk x 19 - 2 

9=0 djiO 



and (hi)ifi respectively. There is automatically no dependence on Qi so the degree zero 
divisor potential is the same as the degree zero free energy (cf. |5Uj): 



F°{t u ... i t k ;x) :=F (t 1 ,...,t k ;x) 

k t 3 t 2 t \ 1 

i=l ' i^j " ijtjjjtl^i / 

k oo 

+ ^(^)i,ot 4 + ^(l) 3 ,o^- 2 . (20) 

i=l 3=2 

Note that degree zero genus 0(1) terms are the only parts of the free energy depending 
on powers of £j rather than just their exponents e*\ When X is compact these terms 
reflect its classical cohomology, namely [50] 



{hihjhi) 0fi = / hiUhjU h t 



x 



(hi) 



1,0 



-1 / hiUaiX). (21) 



In particular, they vanish unless X is a threefold. Higher genus contributions were 
computed in the celebrated paper of Faber-Pandharipande 

(-l)g \B 2g \\B 2g „ 2 \ I 
(20-2)! 2^(2^-2) ' 2 

-l)«- 1 (20-l)S 2gJ B 2fl _ 2 1 



(^O = A- n , of n. ON • n / (<*(*) " Cl (X) U C2(X)) 



(2g-2)(2g)\ 2 



(c 3 (X)- Cl (X)Uc 2 (X)), g>2. (22) 



v 



Here as before are Chern classes and B n are the Bernoulli numbers defined via a 

generating function [5] 



oo „ 



^t^E^ 1 ? ( 23 ) 

— 1 ^— ' n\ 

n=0 

The only non-zero odd-indexed number is Si = —1/2 and B = 1, B 2 = 1/6, B± = 
-1/30, B & = 1/42. 

One sees from (|22p that higher genus contributions all vanish for non-threefolds even 
when non-divisor invariants are taken into account because the Chern classes integrate 
to zero. However, genus 0(1) terms may still survive if X has cohomology classes of 
appropriate degree to cup with c 2 (X) and each other. But the divisor invariants still 
vanish for dimensional reasons. Also note that (|22|) simplifies for Calabi-Yau threefolds 
since C\(X) = and j x c%(X) = is the Euler characteristic of X. Thus, for compact 
Calabi-Yau threefolds 

(20 -2) (2,/)! 2~' 9 - 2 ' (24) 
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When X is non-compact but a^O the moduli M 9 „(X, a) may still be compact. This 
usually happens if geometry forces images of stable maps to stay within a fixed compact 
subset of X, e.g. this is the case for the resolved conifold [26, 33J. Then the virtual class 
is still defined and no new problems arise. However, if a = factorization ffl9|) forces 
M g , n (X, 0) to be non-compact always. To the best of our knowledge no virtual class 
theory exists for non-compact moduli so technically (71 ■ • • 7n) 9 ,o f° r non-compact X are 
not defined at all. 

Leaving the land of rigor and arguing like string theorists we notice that for Calabi- 
Yau threefolds (I24|) still makes sense and can be taken as the 'right' answer even for non- 
compact X. This is consistent with a formal localization computation [8j. Unfortunately, 
for g = 0,1 the invariants contain insertions and we really need to know how to interpret 
the integrals over X in (12"T|) . In physics literature it is suggested that they correspond 
to integrals over "non-compact cycles" [16] that can perhaps be interpreted as duals to 
compact cohomology cocycles [10]. We conclude that for the resolved conifold (x(X) = 2) 
the degree zero free energy has the form 



F° (t- r) - P3 ^ + m (t) + Vf-1 V?' 1 ^ 9 ~ jjjVgggzg r 2g-2 

*0(-i)®0(~i)[ t > x )--^r+Pi[ t ) + 2^ L ) (2g-2){2g)\ ' ^ ' 

where Pi are degree i homogeneous polynomials with rational coefficients. We should 
mention that there are reasonable ways [IB] of assigning values to P3,pi at least for local 
curves (see [11]) from equivariant and mirror symmetry viewpoints. For the resolved 
conifold they yield 

* 0( _i,eo(-i)(*, *) ~ 6 x2 + 12 + 2J, V {2g _ 2 )(2g - 2)! % 

and this function can be recovered from the mirror geometry. However, it appears that 
Donaldson-Thomas and Chern-Simons theories store classical cohomology information 
more crudely. We shall see that in genus 0, 1 this answer or even the general template 
(125]) is inconsistent with exact duality (see discussion after Corollary [1]) . 

Definition 2. The (full) Gromov-Witten free energy is F := F° + F' and the (full) 
Gromov-Witten partition function is Z := exp(F) = Z°Z f , where F',Z f are reduced 
versions from Definition^ As before we write Fx,Zx to indicate the target manifold if 
necessary. 

For the resolved conifold we get from (Tl8|) 

00 

F 0{ _ m0 (-i)(t;x) = -F0(_ 1 ) eO (_ 1 )(t; x) + J>) S)d e dt x 2 °- 2 . (26) 



9=0 
d=l 



The positive degree part converges to a holomorphic function in an appropriate domain 
of t, x (recall that t is a negative Kahler volume). The same holds for all toric Calabi-Yau 
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threefolds and for them the partition function is given directly by the topological vertex 
[21 ED ESI SI] • We will discuss the case of the resolved conifold in more detail in Section El 
But the degree zero part is not so well behaved. The sum in (1251) diverges and fast! By 
a classical estimate for Bernoulli numbers 

^ 9 2 2 9 -i < ™ < ^(22^-1 _!)' KH) 

and the general term in ( 1251) grows factorially for any x^O. Coming up with a space of 
formal power series where the sum lives is neither difficult nor helpful. A helpful insight 
comes from the conjectural duality with Donaldson-Thomas theory [12] that suggests to 
view ( l25"l) as an asymptotic expansion of a holomorphic function at a natural boundary 
point. The function in question is the Mac-Mahon function M(g), the point is q — 1 and 
the relation to (1251) is q = e lx . We inspect this idea in the next section. 



2 Donaldson-Thomas theory 
and the Mac-Mahon function 

In this section we clarify the relationship between degree zero Gromov-Witten invariants 
and the Mac-Mahon function 

oo 

M(g) ^n^-^T"' I?I<1- (28) 

n=l 

This is a classical generating function for the number of plane partitions [1] . [381 1.5.13]. 
More to the point, it appears in [25j 02] in the generationg function of degree zero 
Donaldson-Thomas invariants of Calabi-Yau threefolds. 

Donaldson-Thomas invariants. The Donaldson-Thomas theory provides an alter- 
native to the Gromov-Witten description of holomorphic curves in Kahler manifolds, 
utilizing ideal sheaves instead of stable maps. Intuitively, an ideal sheaf is a collection 
of local holomorphic functions vanishing on a curve. This avoids counting multiple cov- 
ers of the same curve separately and Donaldson-Thomas invariants are integers unlike 
their Gromov-Witten cousins. Counting ideal sheaves is at least formally analogous to 
counting flat connections (i.e. locally constant sheaves) on a real 3-manifold, and the 
Donaldson-Thomas invariants are holomorphic counterparts of the Casson invariant in 
Chern-Simons theory |59j . 

The genus g of a stable map is replaced in a Donaldson-Thomas invariant D K ^ a by the 
holomorphic Euler characteristic k of an ideal sheaf. As conjectured in [42J and proved 
in [35] the degree zero partition function of a Calabi-Yau threefold X is given by 

oo 

Z°x(q) :=J2 D ^q K = M(-q)^ x \ 

K = 
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where as before x(X) * s the classical Euler characteristic. 

Since both kinds of invariants are meant to describe the same geometric objects one 
expects a close relationship between them. Indeed, it is proved in [12] for toric threefolds 
and conjectured for general ones that reduced partition functions of Donaldson-Thomas 
and Gromov-Witten theories are the same under a simple change of variables. This 
equality does not extend directly to degree zero but it is mentioned in [12] that the 
Gromov-Witten F° is the asymptotic expansion of lnM(e i:c )2xW at x = (note the 
extra 1/2 in the exponent). 

A quick look at (1251) tells one that even for the resolved conifold this can be true at best 
for g > 2 since no extra variables are involved in the Donaldson-Thomas function. We will 
see that this is the case but the complete asymptotic expansion involves some interesting 
extra terms that are perplexing from the Gromov-Witten point of view. However, the 
Mac-Mahon factor is exactly reproduced in Chern-Simons theory (Lemma [1]). Moreover, 
with asymptotic expansions one has to specify not just a point but also a direction in the 
complex plane in which the expansion is taken, and the correct direction here is not the 
obvious (real positive) one. 

^-resummation. To avoid imaginary numbers we first consider lnM(e _x ) instead of 
lnM(e lx ). For motivation, we start with a provocative 'computation' that converts an 
expansion in powers of e~ x into one in powers of x for a simpler function: 

n=l n=l k=0 k=0 n=l k=0 

The last two equalities are nonsense of course: the interchange of sums is illegitimate 
and Y^Li = Y^=i nk i s ( ver y) divergent. It certainly does not converge to £(— k) 
for positive k, although by definition £(s) := Y^=i Res > 1 is the Riemann zeta 
function [5]. Nonetheless, the end result is almost correct. Indeed, by definition of 
Bernoulli numbers (T23l) 

C " 1 X 'tB*! (30) 



1 — e~ x x e x — 1 

j=0 

and C(-k) = k > 1; C(0) = -1/2 [5], so (31) 

--Vy^^ = --i-y^^^y(-#/ ( 32 ) 



l-e~ x x 2 ^ ?! x 2 + x ^ v ' k\ 

j=2 J k=l y ' k=0 

In other words our 'computation' (1291) only missed the first term 1/x. 
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A similar feat can be performed with lnM(e x ). First, we compute 

OO OO OO , - nx \k oo -I °° 

In M(0 = - J2 n Ml - e— ) = J>£ }~ = E £ E 

n=l n=l fc=l fc=l n=l 



-kx\n 



f> 1 e' fc * f>l 1 _- csch 2 (f) 

So far all the manipulations are legitimate assuming x > 0, although they would not be 
if we used e lx instead of e~ x . Next recall the Laurent expansion at zero of csch 2 : 

,2/ \ ( 2 # - 1) B 2g 2 q-2 (oa\ 

csch (z) = -2^ tttj, * 9 ■ (34) 

9=0 [A9) - 

One can now pull the same trick as in ([2"§j) of interchanging sums and replacing divergent 
power sums of integers with zeta values. Namely, 



In M(e-) = -Y-Y 4 i2g r ]] B2a 2 2 - 2 ( ^) 



2g-2 



l=l _ v (2g - 1) B 2g v J_ l=l _ v (2g-l)£ 2g C(3-2g) 2 ,_ 2 
^ (2g)\ 2^ k 3-2 9 2^ r 2g )\ 

9=0 v y; fc=l g=0 v y; 

C(3) C(l) , ^ (2^-1)^^-2 ^-2 
-I* 12 + ^ (2, -2) (2g)\ ' (35) 

where we used (1311) in the last equality This series is even more problematic than the 
one in fl29|) which at least made sense and converged for \x\ < 2ir. Now not only does it 
diverge factorially (see fl27|) ) but also £(1) makes no sense at all since ( has a pole at 1. 
Nonetheless, dropping the singular term ^r, the 'infinite constant' — and formally 
replacing x by — ix in the sum, we get exactly the higher genus Gromov-Witten free 
energy in degree zero fl25l) . 

The procedure used in fT29|) . fl35l) can be traced back to Euler and in a more so- 
phisticated guise is used in quantum field theory under the name of ^-resummation or 
^-regularization pQ . The amazing fact is not that this is reasonable to do in physics (one 
can argue that ((—k) has the same operational properties as the non-existent Y2^=i nh )> 
but that it actually produces nearly mathematically correct answers. Unlike a physical 
situation, where a sensible answer is taken as a definition for an otherwise meaningless 
quantity, here we have an identity where both sides make perfect sense (as a holomor- 
phic function and its asymptotic expansion respectively) and only the passage from left 
to right is odious. 

Mellin asymptotics. A fix is a well-known Mellin transform technique that not only 
takes care of singular terms, divergent expansions and infinite constants but even explains 
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Figure 1: Barnes contour for Mellin transforms 



why the double blunder in ( 1291) and ( 135]) computes most of the asymptotic correctly [15] . 
We use it here to make the relationship between the degree zero invariants and the 
MacMahon function precise. Recall that given an integrable function on (0, oo) with a 
possible pole at and polynomial decay at oo its Mellin transform is 



POO 

mtf(s) := / x s - 1 f(x)dx. 
Jo 



The transform is defined and holomorphic in the convergence strip Res G (a, (3), when 
f(x) ~ 0(x~ a ) at and ~ 0{x~P) at oo, assuming a < (3. It is most useful when Wlf 
admits a meromorphic continuation to the entire complex plane since location of the 
poles determines asymptotic behavior of the function at and oo (see [15] and below). 
For example, 97l[e _n:r ](s) := T(s)/n s in Res G (0, oo) extends meromorphically with the 
poles of the gamma function located at s = 0, — 1, —2, .... Analogously, 



971 



1 - e~ x 



r (s 



^aJt[ e -H(s) = EA Z = r ( s K(^) mResG(l,oo) (36) 



n=l n=l 



extends with one additional zeta pole at s = 1. 
The inverse Mellin transform recovers / as 

i pc+ioo 

f( x ) = — -/ Wtf(s)x- S ds for cG {a, /3) (37) 

llli J c—ioo 

assuming absolute integrability along Re s = c. In the cases of interest to us all the poles 
are located on the real axis to the left of a. If the transform satisfies appropriate growth 
estimates one can shift the integration contour in (|37|) to run counterclockwise along the 
real axis from — oo to a and back, Fig. [TJ This reduces (137|) to a sum over residues at 
the poles by the Cauchy residue theorem: 

oo 

f(x) = J2 ^s [3K/(s)]ar* (38) 

n=0 



15 



If 7„ are integers and the series converges / must be real-analytic on (0, oo) with at 
worst a pole at 0, and the residues give its Laurent coefficients at 0. For example, one 
can compute the Taylor expansion of e~ x at using that 9Jt[e _:E ](s) := T(s) and the poles 
— 7„ = — n of T are simple with the residues (— l) n [5]. However, in most cases the series 
( 158]) diverges for all x ^ and (136]) is such a case. Under analytic assumptions that we 
do not reproduce here the following weakening of (|38p is still true 



If —In o^re order m n poles of (meromorphic continuation of) VJlf(s) and its 
Laurent expansions at — j n have the form 

A A m "~ 1 A 

era r / \ ^-nO rt-nl \ ~* ^-nk 



S + ln {S + In) 2 ^ {S + In) 



k+l 



then an asymptotic expansion of f at x = is 

' [nk \n k x)x^. (3!)) 



oo / mn-1 /_-, \k A 

f(x)~Y,[A n0 -A nl \nx+ J2 - J xT 

n=0 \ k=2 



Now it becomes clear where the extra \jx in fl32l) came from. In addition to gamma 
poles in (136]) that produce terms ((—n)x n there is also a simple pole of ((s) with 
residue 1 that gives T(l) ■ 1/x = 1/x. Thus, ( 132]) is at least an asymptotic expansion of 
e~ x /(l — e~ x ) at 0. The fact that it actually converges to the function is a rare bonus. In 
general, even if (139]) does converge it is not necessarily to the original function, see |15j . 
This technique extends to general Fourier sums (or harmonic sums) of the form 

oo 

f( x ) = ^2 a k9(^kx) (40) 

k=l 

because their Mellin transforms can be easily expressed in terms of those of the base 
function g |15j . One can think of them as sums of generalized harmonics with amplitudes 
ak and frequencies Uk, the usual ones corresponding to g(x) = e lx , u>k = ±fc. Indeed, 

°° POO °° POO 

Ttf(s) = J2 ak / x'- 1 g{u k x)dx = ^2-^ / x s - 1 g(x)dx = D(s)Wlg(s), 

k=l J ° k=l U k JO 

where D(s) := YlT=i a k/ u k ls ^ e Dirichlet series of the sum. If D(s) is entire and VJtg(s) 
only has simple poles at s = 0, — 1, — 2,... then 

oo 

f(x) ~ V Res [Wlg(s)]D(-n)x n . 

' J s=—n 
n=0 



If moreover g itself is entire and decays fast enough on R + then g(x) = ^2^=0 9n xn , 

■\oo 
;n=Q • 



Res s= _ n [9Jtg(s)] = g n and f(x) ~ Y^=o9nD{.— n )x n . The same answer can be obtained 
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by a (legitimate under the circumstances) interchange of sums in (j40j) : 



f( X ) = E ak E 9n(u k x) n = 9n I ak W k J % n = E 9n D(-n) l" 

k=l n=0 n=0 \k=l ) n=0 

In particular, this expansion is not just asymptotic but convergent. If D(s) is not entire 
but only meromorphic the last two equalities fail. However, fl39|) still ensures that formal 
interchange of sums gives the regular part of the asymptotic expansion correctly as long 
as D -poles are real-part positive! This is precisely what happened in f[2lj]) . 

Ramanuj an- Wright expansion. The situation in (1351) is more complicated. We com- 
pute from (1551) 

OO ^ OO OO ^ OO -p/ \ 

m [mM(e-)] M = £ i£ »»[r-'1M = £ I £ =2£ 

k=l n=l fc=l n=l ^ ' 

Now assume that Res is large enough for the double series to converge absolutely, e.g. 
Res > 2, and proceed 

00-p/\ oo oo ^ 

= E = roo E • E ftt = ns) as - 1) co> + 1). (41) 

k,n=l 11 = 1 k=l 

The extra zeta poles occur at s — 1, s + 1 = 1, i.e. s = 0, 2 and s = becomes a double 
pole. Formula fl39|) now yields an asymptotic expansion for lnM(e _a: ) that we state as a 
theorem. This is a particular case of asymptotic expansions for analytic series obtained 
by Ramanuj an who used a rough equivalent of Mellin asymptotics, the Euler-Maclaurin 
summation (see [HI Theorem 6.12]). Ramanujan's considerations were heuristic and in any 
case remained unpublished until much later. The first rigorous asymptotic for lnM(e~ x ) 
is due to Wright [63]. We sketch a proof for the convenience of the reader. 

Theorem 1 (Ramanuj an- Wright). Let M(q) := Yl™Al — q n )~ n , \q\ < I be the MacMa- 
hon function. Then lnM(e -x ) has the Mellin transform 271 [In M(e _:r )] (s) = r(s) ((s — 
1) ((s + 1), Res > 2 and its asymptotic expansion at x = along M + is 

In M(e-) ~ ^ + — + C'(-l) + T to- 1 ) 3 ** 3 *?-* x *9-\ (42) 
K 1 x* 12 ^ ; (2g-2)(2g)\ 1 ; 

Proof. Recall that ((s) has "trivial zeros" at negative even integers [5]. Poles of T(s) at 
negative odd integers are therefore canceled by zeros of ((s — 1). Analytical assumptions 
needed for (1391) to hold are satified here by the classical estimates for F and ( [5j. The 
contributing poles are: 

(i) Gamma poles at s = —2, —4, . . . , —2g, . . . with residues ^$J- £(—2g — 1)C(1 — 2p); 
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(ii) Simple pole of ((s — 1) at s = 2 with residue 1 • r(2)£(3) = C(3)i 

(iii) Double pole of T(s), ((s + 1) at s = 0. 
We have from the first two items and (1311) 

f + | ^CM, -!){(!- 2,) ^ = f + g \!!^-' 

To take care of the double pole we need more than just the residue. By the well-known 
properties of T and £ 



2^ 



r(l + s) = 1 - 7s + 0(s 

C00 = -^-+7+0(3-1), 
s — 1 

where 7 is the Euler constant. Thus 

r(a) C(a + 1) = + + 1} = - 7 + OOO) (~ s + 7 + O(a)) = ^ + 0(1), and 

i» c(« - 1) C(s + 1) = (± + o(i)) (c(-i) + C'(-i> + o(, 2 )) 

S =C^ + ^ + (i) = =l^ + ^ + (i). 

s z s s z s 

By (1391) the corresponding terms in the asymptotic expansion are lnx/12 + ('(—1) and 
it remains to combine the expressions. □ 

In hindsight, it is amusing how much of (H2|) is visible in the naive expression (|35|) : 
not just the regular part but also ((3)/x 2 and even 1/12 in front of the logarithm. The 
only hidden term is ('(—1), sometimes called the Kinkelin constant [I], and for this reason 
perhaps it is usually missing in physical papers. 

Stokes phenomenon and the natural boundary. As already mentioned the rela- 
tionship between q and x is q = e tx not q = e~ x . Replacing formally x by — ix in (H2l we 
recover the infinite sum of (|25[) along with three extra terms 

C(3) \n(-ix) 
x 2 12 ^ l h 

How legitimate is this substitution? Had (H2l been a convergent Laurent expansion there 
would be no such question. But it is asymptotic and represents lnM(e~ x ) only up to 
exponentially small terms (more precisely, "faster than polynomially small" but we follow 
the standard abuse of terminology). It is well-known that such expansions depend on a 
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direction in the complex plane in which they are taken. As one crosses certain Stokes 
lines originating from the center of expansion exponentially small terms may become 
dominant and change the expansion drastically. This change is commonly known as the 
Stokes phenomenon. Moreover, for an asymptotic expansion in some direction to exist 
the function must be holomorphic in a punctured local sector containing this direction in 
its interior. Switching from x to — ix while keeping x real positive forces us to approach 
q = e l ° = 1 along the upper arc of the unit circle, i.e. along a purely imaginary direction. 
For an asymptotic expansion in this direction we need to have M(g) analytically continued 
beyond the unit disk \q\ < 1. But can it be continued? 

The definition (1281) does not look very promising. In fact, it strongly suggests that 
M(q) has a singularity at each root of unity. But roots of unity are dense on the circle 
making it a natural boundary for M(q) and no analytic continuation exists. It turns out 
to be quite hard to turn this observation into a proof but Almkvist shows [I] that if a/6 
is a proper irreducible fraction then 

lnM^V*) ~ M + Ai nx + o(l) 
b 6 x z 12 

for real positive x. Thus, every root of unity is indeed singular and \q\ = 1 is the natural 
boundary. 

This forces us to reconsider keeping x real in lnM(e ja: ). Should x approach from 
the positive imaginary direction we can set x = iy with y > and Theorem [1] gives us an 
asymptotic expansion in y. We can rewrite it as an expansion in x of course as long as 
it is understood that x in it is positive imaginary. This may seem like an underhanded 
trick but it is not. The natural domain of lnM(e ja: ) is the upper half-plane and the only 
distinguished direction in its interior is the positive imaginary one. 

Corollary 1. Asymptotic expansion of lnM(e llE ) at x = along i~R + is (taking the 
principal branch of the logarithm) 



lnM(e i: 



C(3) ln(-ix) lU f ( lVq _ 1 (2ff- l)B 2g B 2g _ 2 2 2 

^ + ^2- + C( - 1) + y- 1} (2g -2) (2,)! X 

C(3) hix _ tu ^ 1 (2g - 1) B 2g B 2g . 2 2c/ _ 2 

x 2 + 12 +U ' 24 + ' {2g-2){2g)\ ' 1 > 



Comparing this to (125]) one ought to be somewhat perplexed. If we are to take fj43l 
at face value then ps(t) = — C(3), Pi(t) = C'(~ 1) — 7ri/24 (?!) and there is no space for 
lnx/12 at all. Aside from the fact that pi-s are supposed to be homogeneous polynomials 
of the corresponding degree the numbers involved are not even rational, £(3) by Apery's 
famous result. Nevertheless, the MacMahon factor appears as is in the Chern-Simons 
partition function, see Lemma [TJ 

The disappearence of extra variables and appearence of irrationals suggests that some 
kind of averaging is involved. It would not explain In a?/ 12 but we may guess that aver- 
aging of pi(t) is divergent and has to be regularized giving rise to an anomalous term. 
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Why Donaldson-Thomas theory does not reproduce the degree zero contributions in low 
genus is beyond our expertise. However, from the Chern-Simons vantage point this ought 
to be expected. The idea of large N duality is that the same string theory is realized 
on manifolds with different topology [HI HE]- However, the degree zero terms in genus 
0, 1 are exactly the ones that record the classical cohomology of the target manifold, 
see f l2~Tj) . Although some relation between topologies of manifolds supporting equivalent 
string theories may be expected, the entire cohomology ring is certainly too much to sur- 
vive a geometric transition. Therefore, these classical terms can not enter an invariant 
partition function except via averages that remain unchanged by such transitions. 

3 Topological vertex and 

partition function of the resolved conifold 

This section and the next are to be read in conjunction. We review the salient points 
of two combinatorial models, the topological vertex [2, ED ESI HI] and the Reshetikhin- 
Turaev calculus [HI EE], highlighting the differences but more importantly the parallels 
between them. The former computes the Gromov-Witten invariants of toric Calabi- 
Yau threefolds and the latter computes the Chern-Simons invariants of all closed 3- 
manifolds. The reason to compare them is the conjectural large N duality between the 
two. Both models encode their spaces into labeled diagrams and then assign values to 
them according to Feynman-like rules. However, the encoding and the rules are quite 
different despite intriguing correspondences. The reason we use the topological vertex 
instead of just summing up (1261) as in [H] is that it directly gives the partition function in 
correct variables and in an appealing form. Comparing the answer to the Chern-Simons 
one it becomes reasonable to express it in a closed form via the quantum Barnes function 
(Theorem E]). 

Toric webs. Just as the Reshetikhin-Turaev calculus [8j [58] the topological vertex is 
a diagrammatic state-sum model. This means that geometry of a space is encoded into a 
diagram, a graph enhanced by additional data, and the value of an invariant is computed 
by summing over all prescribed labelings of the diagram. In the Reshetikhin-Turaev 
calculus the diagrams are link diagrams representing 3-manifolds via surgery [521 [58]. In 
the topological vertex they are toric webs representing toric Calabi-Yau threefolds. 

A toric web is an embedding of a trivalent planar graph with compact and non- 
compact edges into K 2 that satisfies some integrality conditions [211 ESI HI]- Namely, 
vertices have integer coordinates and direction vectors of edges can be chosen to have 
integer coordinates. Moreover, if the direction vectors are chosen primitive (without a 
common factor in coordinates) any pair of them meeting at a vertex forms a basis of 
Z 2 and every triple at a vertex if directed away from it, adds up to zero. Examples for 
the resolved conifold 0{-l) © O(-l) and the local CP 1 x CP 1 (i.e. the total space of 
T 1,0 (CP 1 x CP 1 )) are shown on Fig. [2J where the primitive directions of non-compact 
edges are also indicated. Toric webs related by a GL-{L transformation and an integral 
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Figure 2: Toric webs of 0{-l) © 0(-l) and local CP 1 x CP 




Figure 3: Toric graphs of 0(-l) © £>(-!) and local CP x CP 



shift represent isomorphic threefolds. For this reason we did not label the vertices on 
Fig. |2l one may assume that one of them is (0, 0) and all compact edges have the unit 
length. The toric web is a complete invariant of a toric Calabi-Yau. Indeed, the moment 
polytope of the torus action can be recovered from it [361 4.1] and therefore the threefold 
itself up to isomorphism by the Delzant classification theorem [55J. Analogously, a 3- 
manifold is recovered from its link diagram up to diffeomorphism by surgery on the link 
[52l [58] . Having toric webs rigidly embedded in R 2 is inconvenient; one would prefer to 
treat them as abstract graphs, perhaps with additional data. This is possible at least as 
far as the topological vertex is concerned although the resulting graphs may no longer 
be complete invariants. 

Tracing back the construction of a threefold from its web one concludes that the 
vertices correspond to fixed points of the torus action and compact edges correspond to 
fixed rational curves (copies of CP 1 ). Being rational curves sitting inside a Calabi-Yau 
threefold their normal bundles are isomorphic to 0(n — 1) © 0(—n — 1), n = ±1, ±2, .... 
The framing number n e for each edge e is assigned the value from the normal bundle 
type of the corresponding curve. This only determines n e up to sign and the edge must 
be oriented to specify it. Although on their own these orientations are chosen arbitrarily 
they must be aligned with the framing numbers, the exact rule is given in [36j 4.2]. 

If £i, . . . , is an integral basis in H 2 (X, Z) as in Section [1] then each edge curve C e 
represents a homology class expressible as a linear combination [C e ] = ?7ii£i + • • • + m^k, 
rrii G Z. One requires these homology relations to be attached to the edges as well. The 
result is a graph called the toric graph. Toric graphs for 0(—l) © 0(—l) and the local 
CP 1 x CP 1 are shown on Fig. [3l Framing numbers and homology relations are the only 
data aside from the topology of the web used in the topological vertex. We emphasize 



21 



X= (4,2,2,1,0, ... ) 



K= (4,3,1,1,0, ... ) 



Figure 4: Young diagram and its conjugate 

that both can be recovered algorithmically from the web itself without any recourse to 
the original threefold [23], [361 4.1]. 

Partitions and Schur functions. We wish to briefly describe the topological vertex 
algorithm to see how g-bifactorials naturally emerge from it. This requires some basic 
information about partitions [35] that appear in the Reshetikhin-Turaev calculus as well. 
Partitons serve as labels in state sums defining the invariants. A partition A is an element 
of with only finitely many non-zero entries that are nonincreasing, i.e. 

A = (Ai, A 2 , . . . , X N , 0, . . . ), A i eZ + , X x > X 2 > ■ ■ ■ > X N > 0. 

Let V denote the set of all partitions. The number of non-zero entries /(A) is called 
the length of a partition and the sum of all entries |A| := Ai + A2 + ■ ■ • + Aat is called 
its size (or weight). Partitions are visualized by Young diagrams, rows of boxes stacked 
top down with Aj boxes in i-th row, Fig. HI The conjugate partition X' is obtained 
visually by transposing the Young diagram along the main diagonal and analytically as 
A^ := max{j|Aj > i}. Note that A" = A and X[ = 1{X), \X'\ = \X\. Another relevant 
characteristic of a partition, sometimes called its quadratic Casimir, is 

00 

H (X) :=5^Ai(Ai-2i + l), x(A') = -x(A), x(A) G 2Z. (44) 
i=i 

Partitions represent possible states of compact edges in a toric graph and a combination of 
partition labels for each edge represents a state of the graph [31]. The partition function 
is then obtained by summing over all possible states. 

Amplitudes (see Definition HJ) of a labeled graph are defined via a specialization of 
Schur functions s\ indexed by partitions. They are symmetric 'functions' in the sense 
of Macdonald [38], i.e. formal infinite sums of monomials in countably many variables 
that become symmetric polynomials if all but finitely many variables are set equal to 
zero (more technically, if monomials containing any variable outside of a finite set are 
discarded from the sum). For instance, if A = (l n ) := (1, 1, . . . , 1, 0, . . . ) then s^ is the 

n times 

n-th elementary symmetric function: 

s {1 n)(x) = e n (x) := ^2 x h--- x i n - (45) 

l<ii<— <i n <oo 
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In general, s\ are polynomials in the elementary symmetric functions given by the Jacobi- 
Trudy formula s\ = det(e A ^_j + j), 1 < i,j ' < /(A') = \\. For example, 



'(2,1,0,...) 



X 



oo 

e2 e ' 



e 3 ei 



eie 2 - e e 3 = 22xi- 22 XiXj - ^ XiXjX k . 

i=l i<j=l i<j<k=l 



Since e n are homogeneous of degree n the Jacobi- Trudy formula implies that s\ are also 
homogeneous of degree |A|, i.e. s\(ax) = a^ x ^s\(x). Moreover, s\, A G V form a linear 
basis in the space of symmetric functions, in particular s\s^ = J2vev c uv s v It turns out 
that are non- negative integers that vanish unless \v\ = |A| + They are the famous 
Littlewood-Richardson coefficients [38J. 

Specializations of Schur functions appearing in the topological vertex are obtained by 
specializing the formal variables Xi to elements of a geometric series possibly modified at 
finitely many entries. Such specializations were extensively studied by Zhou [M]- Define 
the Weyl vector p by 

13 , / . T " 



Note that p is not a partition. Introduce a new formal variable q and for any vector £ 
set g£ := (g^ 1 , g^ 2 , . . . ) so in particular g p = (g~2 ( g~2 ( . . . ) is a geometric series. 



Definition 3. One-, two- and three-point functions of the topological vertex are respec- 
tively JM [M 

Wx(q) :=s x (q p )=Wxo(q) 

VM?) ■■= s x (q P K(q X+p ) (46) 

There is a shorter expression for the three-point function via the skew Schur functions 
[3Tj 161] but we do not need it here and (|46p is somewhat reminiscent of the Verlinde 
formula in Chern-Simons theory [61]. We assume q G C\R_ and qz is then defined 
by the principal branch of the square root. One can see by inspection from (|45|) that 
e n (<? A+p ) converges for \q\ > 1. Since the Schur 'functions' s M are polynomials in e n they 
are also well-defined as honest functions of q upon specializing to q x+p . 

To be consistent with the usual basic hypergeometric notation [T5] we wish to switch 
from \q\ > 1 to |g| < 1. This can be done using a symmetry of the two-point functions 



mM = (-i) |a|+m Wav(? _1 ) = (-i) lM+M sy(q- p )sAq~ x '~ p )- (47) 

This identity is a curious one since the two sides never converge simultaneously (both 
diverge for \q\ = 1). It has the same meaning as a more familiar identity 



oo oo 



i=l * * i=0 i=l 
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—* 

Figure 5: Partition triple at a vertex A„ := (A',//, v) 



where the two sides never converge simultaneously either. In fact, Wx^q) are rational 
functions of and can be analytically continued to C\R_, ( I4T1) expresses this analytic 
continuation. 

The appearence of g-bifactorials in partition functions is due to the Cauchy identity 
for Schur functions [38j 164"] 

oo 

Y, s ^ x ) s ^y) uW = n^ + ux lVj ). (48) 
xev i,j=i 

If Xi = q 1 " 1 , yj = qi~ l the righthand side of (|4"8]) becomes 

oo oo 

n (i +« 9 i -y- 1 ) = n (! + u <i t+j ) = («; ( 49 ) 

ij=l j,j=0 

Note that although (HHj) is a formal identity if both sides converge as in (H9|) it holds as 
a function identity. 

Partition functions as state sums. Let us now inspect the state sums appearing 
in the topologial vertex. Let V and E c denote the sets of vertices and compact edges 
of a toric graph respectively. Choose an arbitrary orientation for each element of E c , 
this determines the sign of the framing numbers. Assign a formal variable a« to each 
element of a basis ^ G H 2 (X, Z) and set a e := a™ 1 • • ■ a™ k for the corresponding edge 
curve [C e ] = mi^i + ■ • • + Finally label all compact edges by arbitrarily chosen 

partitions A e G V and non-compact ones by the trivial partition G V . Triples of 
partitions A„ := (Ai,A 2 ,A 3 ) are then assigned to each vertex according to the following 
rule 

Starting with any of the three edges and going counterclockwise around the 
vertex pick the edge's label if the arrow on the edge is outgoing and its 
conjugate if the arrow is incoming, Fig. [5] 

Non-compact edges present no problem since 0' = 0. This determines the triple up to 
cyclic permutation which is enough since := Wa 1 a 2 a 3 has cyclic symmetry. 
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Definition 4. Amplitude of a labeled toric graph relative to a basis £1, . . . , G H 2 (X, Z) 
is given by J2H ESf 

|Ae|(n e +l) a |A e | _ 

The main result of [36! can be stated as follows 



A {Xe} ( ai , ...,a k ;q):=H (-1)1^1(^+1) q ^ a M . J] y^ (g) . (50 ) 



Theorem 2. T/ie reduced Gromov-Witten partition function of a toric Calabi-Yau three- 
fold X relative to a basis £i, • • • , € H2(X, Z) is <n'i> en fry a state sum 

Z x (a 1 ,...,a k ;q) = ^2 ^{A e }(oi, • • • , a k ; q) = A Xl ,...,x lEc{ {a%, . . . , a k ; q) 

i x e}\ ee v Ai,...,A| Bc |e'P 

(51) 

assuming in the second sum that the edges are numbered and Aj := X £i . 

Partition function of the resolved conifold. Here H2(X, Z) is one-dimensional and 
f = [CP 1 ]. There is only one a variable and only one edge. The amplitude A\ for A G V 
is (see Fig. [3] and flS}) 

A\(a; q) : = (-1)I A I ■ W xoo (q) Wx>oo(q) = (-a^W^q) Wx>(q) 

= (_ a )|A| ( _ 1) |A| + |A'| SA/(g -p )sA(g -p ) 

Recalling that |A'| = |A|, — p = (i — |)°!\, and s\ is homogeneous of degree |A| we compute 
further 

= (-apsy(q^)s x (q^) = (-a^q^ s x ,((q*))s x ((q>)) = (-og-^VCfo*))^)). 

Suppose that a is small enough for ^2 Xe -p A x (a; q) to converge then we get by Theorem [2] 
and the Cauchy identity f|48l) 



Z' x {a;q) = ™M Wv(g) (~a) W = ^ s x ,((q%s x ((q>))(-aq-T (52) 

CO CO 

= n (i + m = n o- - ^<i i+3 ) = im ?)s (53) 

i,j=l i,j=0 

If we accept the MacMahon function as the degree zero partition function of the resolved 
conifold (despite the issues discussed after Corollary [1]) then 

Z x (q)=M(q) Xj f ' 1 



ir=i(l-^ (q;q)£ y 
We conclude that the (full) Gromov-Witten partition function of the resolved conifold is 

(gq;q) c 2) 
(q;q)L 

as used in the Introduction. 



Z x {a- q) = Z° x (q) Z' x (a; q) = (54) 
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4 Reshetikhin-Turaev calculus and 
partition function of the 3-sphere 

As explained in the beginning of the previous section this one is complementary to it. We 
briefly review the sl^C Reshetikhin-Turaev calculus P, 58] in a form that invites analogies 
with the topological vertex. In particular, we forgo the usual terminology of dominant 
weights and irreducible representations of sl^C and rephrase everything directly in terms 
of partitions. The immediate goal is to compute the partition function of S 3 in a suitable 
form and compare it to the one for the resolved conifold (Theorem [3]). 

Whereas computation of Gromov-Witten invariants in all degrees and genera is an 
open problem (beyond the cases of toric Calabi-Yau threefolds [35] and local curves 
[TT]). Reshetikhin-Turaev calculus provides an algorithm for computing Chern-Simons 
invariants for arbitrary closed 3-manifolds, especially effective for Seifert-fibered ones 
[10] . This circumstance combined with explicit large N dualities for toric Calabi-Yau 
threefolds is the secret behind physical derivation of the topological vertex. To be sure, 
there is a catch. The Reshetikhin-Turaev model (or equivalently Atiyah-Turaev-Witten 
TQFT [?J [58]) is not a single model but a countable collection of them, one for each 
pair of positive integers k, N known as level and rank. This would not be much of a 
hindrance if not for the tenuous connection between invariants for different k and N. As 
a rule, geometers are interested in asymptotic behavior for large k [53J, and physicists 
in both large k and large N behavior. Reshetikhin-Turaev sums with ranges depending 
explicitly on k and N are not exactly custom-made for those types of questions. In fact, 
they require significant work even in simplest cases to be converted into asymptotic- 
friendly form. No general method exists; most common ad hoc procedures use Poisson 
resummation [3D] or finite group characters [8| [29]. 

The idea of the Reshetikhin-Turaev construction (related but different from Witten's 
original one [61J as formalized by Atiyah [7J) is to combine some deep topological results 
of Likorish- Wallace and Kirby with the representation theory of quantum groups [8l 158]. 
A theorem of Likorish and Wallace asserts that any closed 3-manifold can be obtained 
by surgery on a framed link in S 3 [52]. This is complemented by Kirby's characteriza- 
tion [27J of links that produce diffeomorphic manifolds as those related by a sequence of 
Kirby moves: blow up/down and handle-slide. Blow up/down adds/removes an unknot- 
ted unlinked component with a single twist and handle-slide pulls any component over 
any other one, Fig. [61 Thus, if one can find an invariant of framed links that remains un- 
changed under Kirby moves it automatically becomes an invariant of closed 3-manifolds 
via surgery. 

Hopf and twist matrices. Framed links can be represented up to isotopy by plane 
diagrams with under- and over- crossings and twists as in Fig. [6] providing a combinatorial 
model of 3-manifolds. Slicing a link diagram bottom to top and avoiding slicing through 
cups, caps, twists or crossings one gets arrays of basic elements Fig. [7J stacked on top of 
each other. 

This decomposition fits nicely with structure of a linear representation category: plac- 
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Figure 6: Blow up/down and handle slide over a trefoil knot 



Figure 7: Basic elements and slicing of a Hopf link 

ing elements next to each other corresponds to tensoring and stacking corresponds to 
composition. It remains to find an object with representation category meeting all the 
invariance requirements. It turns out that it is extremely hard to find one producing 
non-trivial invariants. Classical Lie groups and algebras do not work unfortunately. One 
has to deform the universal enveloping algebras of say sInC into quantum groups and 
then specialize the deformation parameter q to a root of unity q = e 2m /( k + N ), As if that 
were not enough the tensor product of representations has to be modified as well. The 
end result [8j [29J, [58] is a representation-like category with only a finite number of irre- 
ducible representations. For sl^C at level k they are indexed by partitions with Young 
diagrams in the (N — 1) x k rectangle, i.e. 

V k N _ x := {A e V | /(A) < N - 1, Z(A') = Ai < k}. 

In the equivalent language of dominant weights this corresponds to the weights in the 
Weyl alcove of the Cartan-Stiefel diagram of sIatC scaled by k, see [8]. The Reshetikhin- 
Turaev invariants are computed as state sums over labelings of a link diagram with each 
link component labeled by a partition from V^_ 1 , a finite set. 

Thus, unlike in the topological vertex where sums are taken over all partitions and 
are infinite, in the Reshetikhin-Turaev calculus sums are finite with explicit dependence 
on k, N. Once a diagram is labeled morphisms between irreducible representations and 
their tensor products are assigned to the elements from Fig. and then assembled by 
tensoring, composing and eventually taking traces (corresponding to caps) to obtain 
numerical invariants. The hardest ones to compute are the crossing morphisms for they 
depend on the so-called R- matrix of a quantum group [3J [58] . Good news is that for a large 
class of 3-manifolds, the Seifert-fibered ones and others, the use of crossing morphisms 
can be avoided entirely in computing the invariants [Id El] (but not in proving their 
invariance). In terms of conformal field theory they are completely determined by fusion 
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rules without involving the braiding matrices [2^]. This means that the only algebraic 
inputs are the Hopf and twist matrices S and T: 

Sx, = S 00 W Xfl ; T v = T 00 V ■ (55) 

The notation is as follows: 

Wa m is the normalized quantum invariant of the Hopf link Fig. [7J with components 
labeled by partitions A,/i (see more below); 

fi* is the partition sl^-dual to fi, /i* := \i\ — pN-i+i for 1 < i < N — 1 and \x\ := 
for i > N (not to be confused with the conjugate partition /i'); 

X N denotes the gl N coordinates of (N — 1) x k partition A, A^ := \f — | A | / TV, in 
particular pf := \{N - 2i + 1); 

^(A^) := A^ • (\ N + 2p N ) is the quadratic Casimir of A as a dominant weight 
closely related to the quadratic Casimir x(A) of a partition A: 
C 2 {\ N ) = x(A) + iV|A| - |A| 2 /iV. 

We shall say more about the normalization constants Soo, T o below. 

Formulas for were originally obtained by Kac and Peterson in 1984 in the context 
of affine Lie algebras. Their relevance to the Chern-Simons theory was discovered by 
Witten [61]. In 2001 Lukac [37] realized that they are specializations of Schur functions 
of iV variables, namely 

W Xll (N;q) := s x (q pN ) St ,(q xN+pN ) with q = e*& . (56) 

This should be compared to the two-point functions (HS"]) of the topological vertex. This 
realization led among other things to the physical derivation of the topological vertex, 
where W\^(q) are obtained as some loosely interpreted limits of W\^(N;q) [3l [23l WL\ . 
Let us emphasize the differences though: in fj4"6"|) q is a formal variable whereas in (1561) it 
is a number. Moreover, the number of variables in s\, s M before specialization is infinite 
whereas here it is iV < oo. This last circumstance dramatically simplifies many com- 
putations with Wa^ compared to their analogs with because infinite specializations 
often have nice analytic expressions |38j . 

The twist matrix T also has a vertex counterpart in the form of the framing numbers 
n e that contribute factors of q n ^ x ^)l 2 to the amplitude ( |50l) . Incidentally, this explains 
their name. In the Reshetikhin-Turaev sums factors account for twists in link di- 
agrams that in their turn represent framing of a link, i.e. trivialization of its normal 
bundle up to homotopy. If one thinks of strands as thin ribbons the signed number of 
twists gives exactly the signed number of full twists in a ribbon. 

The two-point functions are symmetric = and the one-point functions W\ : = 
W\o = Wox are called the quantum dimensions of representations indexed by A. The 
quantum diameter is 

V 2 := W ^ N ^f (57) 
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and S 00 is simply its inverse S 0Q := V^ 1 . Analogously, T 00 is the inverse of the charge 
factor £ := e 27rlC//24 , where c := k dimc(si^C) / (k + N) is the so-called central charge 
from conformal field theory [29]. Thus explicitly, Too := C _1 = q~ k( - N _1 )/ 24 . These 
normalizations are needed to make the S and T matrices satisfy the defining relations of 
SL 2 Z 

(ST) 3 = S 2 , S 2 T = TS 2 , S 4 = I. 

Note that unlike and q\ C2 ^ xN ^ that depend on k only via q = e 2m ^ k+N ^ this is no 
longer the case for the normalizing constants Sqo, T 00 and this causes major problems in 
relating Chern-Simons expressions to Gromov-Witten ones [8]. 

Reshetikhin-Turaev invariants. Let J\ u ...,\ n {L) denote the colored HOMFLY poly- 
nomial of an n-component link L, i.e. the amplitude of the link diagram computed as 
outlined above after labeling the link components by partitions (colors) Ai, . . . , X n , see 
[HI EH] for specifics. Then the Reshetikhin- Turaev invariant of the 3-manifold M surgered 
on L from S 3 is 

t(M) : =C*"Dr«- 1 J2 Jx 1 ,...,x n (L)W Xl ---W Xn (58) 

Ai,...,A n GP* r _ 1 

where a is the signature of the linking matrix of L [52] . Since S 3 can be obtained from 
itself by surgery on the empty link with J(0) = 1, n = components and linking 
matrix we have 

Z s3 := r(S 3 ) = V- 1 = S 00 = ( ^A(iV;g) 2 ) (59) 

This is the Chern-Simons partition function of the 3-sphere to be identified after Witten 
[62] and Gopakumar-Vafa [21J with the string partition function of T*S 3 . Note that 
( I57p bears some resemblance to the expression (1521) for the reduced partition function of 
the resolved conifold especially if we interpret (— a)' A ' as a convergence factor needed to 
extend the sum from V%_ x to all partitions. However, unlike (1521) that gives Z' x directly 
formula f[5Tj) gives Z~ 2 , so naively one does not expect these partition functions to be 
nearly equal. 

The normalization adopted in (158!) is the one that makes r(M) into an honest in- 
variant. It is due to Reshetikhin and Turaev and is known to differ from the physical 
normalization of Witten [21] coming from the path integral. Although the physical 
normalization gives better agreement with Gromov-Witten theory [32] it has not been 
consistently defined in general. In the known examples it is discerned by comparing t(M) 
to heuristic path integral expansions [53J. 

Partition function of S 3 . Although we used q above as much as possible one should 
not forget that in the context of Chern-Simons theory this is simply a shorthand for 
e 2m/{k+N) anc j a jj relevant quantities are only defined for positive integers k, N. Large 
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N duality predicts among other things that it should be possible to interpret r(M) as 
restrictions of holomorphic in q functions to these special values, perhaps up to some 
parasite factors steming form misnormalizations. It is these holomorphic functions that 
should correspond to the Gromov-Witten generating functions in q. According to this 
philosophy we should try to transform the righthand side of floTj) into an explicit function 
of q and iV so far as possible. A clean way of doing this is not known to the author, 
the problem being the rogue range of summation V^_ 1 . The known ways via Poisson 
resummation [40j 4.2] or group characters [8j Theorem 10] are rather down and dirty and 
we do not reproduce them here. An intermediate answer is 

V^ = (-l) m ^ R Nik + N)"' 1 H [qk^-q^-^y 2 (60) 

l<i<j<N 

and it fulfils our wish only partially. One could replace k + iV by \n.q/2m but this does 
not lead to anything useful. 

In fact, N(k + iV)^ -1 is the volume of the fundamental parallelepiped of a rescaled 
root lattice A r of sl^C, namely of (k + iV) 1//2 A r . Ooguri and Vafa [49, 2.1] replace it with 
the volume of Un in their normalization. Thus, we ignore N(k + iV)^ -1 as a parasite 
factor and focus on the product following it instead. Since J2i<i<j<N^ = N(N — l)/2 
and £i<«,<*0' - i) = N(N 2 - 1)/12 we have 

II (qtM - q-tM) = q- 2 *^ ]J (l - q j ~ i ) . 

l<i<j<N l<i<j<N 

Note that q 21 = ^- 2N / k an( j a power of ( appears explicitly in the Reshetikhin- 
Turaev normalization (1581) . We get for the partition function 

„ N(N—1) N(N 2 -1) , , r 1 . _I 1— T . . 

Z s3 := t(S 3 ) = t^-^q-^r^ [N{k + N) ) 2 ] [ (l - q>~ 1 ) (61) 

l<i<j<N 

The first two factors are 1 in absolute value and can be regarded as framing corrections 
and we already discussed the third (volume) factor. In any case, it turns out that only 
the last product is relevant in the context of large N duality as we now demonstrate 
(cf.|23l Appendix]) 

Lemma 1. For any \q\ < 1 

w-i / AT+l \(2) 

Hi 1 ' ^) = II C 1 - ^ = 1) . (62) 

i<i<j<N n=i \q->q)oo 

where (a;g)oo := (a; 5)00 is the usual q- shifted factorial J7E/ and (a;q)oa was defined in 
the Introduction (jSj). 
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Proof. First we arrange the factors according to the powers of q 



N-l N-l 



Yl (1 - q^) = Yl (1 - q n )^ hi\™+i<N} = JJ (1 - g")^"" , 

l<i<j<N n=l n=l 

where 1$ denotes the characteristic function of a set S. A finite product can be written 
as a ratio of two infinite ones as long as the latter converge, in particular 



TT fl _ a n )N -n = UZl (1 ~ <ffZ = JElMf"" 



□ 

Setting a = q N in f[6"iZj) and comparing to the partition function of the resolved conifold 
( 154)) we see that Chern-Simons theory reproduces both the reduced partition function 
and the Mac-Mahon factor. It also produces an additional one (q; q)^ aside from the 
parasite factors discussed above. This extra factor appears naturally though in the 
quantum Barnes function G q . Its most characteristic property is the functional equation 
G q (z + 1) = Tq(z)G q (z), where T q is the quantum gamma function of Jackson [5]. T q is 
a g-deformation of the classical Euler's T and the classical Barnes function satisfies the 
same equation with g-s removed [45] . As we derive in the next section (Theorem H|) to 
have this property G q must be the product 

G q (z + 1) = (1 - g)-V* (g; q)^ 19 ' , 

i.e. reproduce the righthand side of fl62l) with z = A up to a power of 1 — q. Emergence 
of the same expressions from a simple functional equation is quite intriguing as well as 
the fact that they are the ones common in Chern-Simons and Gromov-Witten theories. 
We summarize the observations made in the last two sections as a theorem. 

Theorem 3. Quantum Barnes function G q is the common factor of the partition func- 
tions of the resolved conifold X and the 3-sphere, namely 

(1 - q)^ Z x (q z ; q) = (q; q)~* G q (z + 1); (63) 



z-l 

2 



1-g )— Z s z (z;q ) = i—q ~ z 2 I -1 \ G q (z + 1). 64 

\2m J 

In (1631) g is the usual variable of the topological vertex and a := q z is the degree variable 
J3b^ . In ([Mil g = e 2m ^ k+N ^ and z = N, where k,N are level and rank. Equivalently, if t 
is the Kahler parameter and x the string coupling constant of JJ^j then q = e lx , q z = e~ l . 
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5 Quantum multigamma hierarchy 



In this section we give a self-contained introduction into the theory of the quantum 
Barnes function and its higher order analogs, quantum multigammas. The main result 
is the alternating formula (1761) for G^f 1 that displays its graded product structure. Along 
the way we describe connections to other special functions that came into the spotlight 
lately, g-shifted multifactorials and quantum polylogarithms. 

Quantum multigammas. Quantum multigammas emerge naturally if one iterates 
two classical constructions. One is the construction of factorials from natural numbers 
with subsequent analytic continuation to complex values via the functional equation 
r(z+l) = zT(z). The other is a g-deformation of T first performed by Jackson, thorougly 
forgotten and then revived by Askey in 1970-s, see [5j. Euler's construction was iterated 
by Kinkelin in 1860-s who turned T(z) into a new z, i.e considered the equation G(z + 
1) = T(z)G(z) but only for positive integers z. Barnes in 1900-s introduced G(z) from 
entirely different considerations that define it for complex values directly and lead to a 
hierarchy of functions with G^(z) = z,G^(z) = T(z),G^(z) = G(z), . . . ; G^(z + 1) = 
G^ _1 )(z)G^(4 Nowadays G(z) is known as the (classical) Barnes function and G^ as 
multigamma functions jUJ . 

It is clear that the hierarchy of functional equations together with a normalization 
= 1 uniquely define G^ on all natural numbers. Extension to complex values can 
not be unique even for V because there are entire functions that vanish on all integers, 
sin(7T2;) for example. However, Bohr and Mollerup proved in 1920-s that if a log-concavity 
condition is added 4^T{x) > for x > then the Euler's T is the only possibility [5]. 
In 1963 Dufresnoy and Pisot generalized the Bohr-Mollerup existence and uniqueness 
result to a wide class of functional equations of the type f{z + 1) — f[z) = <f>(z). It 
constructs not only G^ d > but also the Jackson's deformation of Euler's Y that satifies 
T q (z + 1) = (z) q T g (z) with (z) q := called the quantum number [181 EE]- For higher 
order multigammas the log-concavity condition has to be replaced with log-positivity of 
order d+1, i.e. lnT(x) > for x > 0. After g-deformations of Barnes' multigammas 
appeared in the context of integrable hierarchies [T7j Nishizawa realized that the same 
iteration works for them as well [45] . In other words, there exists a unique hierarchy of 
meromorphic functions satisfying 

(i) GjO(l) = l, Gf\z)= 1 -^ 

H H 1 — q 

(ii) Gf (z + 1) = G^iz) G q d \z) (65) 

d d+1 

(iii) lnG v( x ) ^ for a; > 0, < g < 1 

The last condition is only required to hold for real positive q but it is understood that 
Gq are continued analytically to other values of q. 
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Definition 5 (Quantum multigammas) . The functions G q defined by (165]) are called q- 
multigamma functions, in particular V q := G q is the Jackson's quantum gamma function 

(2) 

and G q := G q is the quantum Barnes (or q-Barnes) function. 

For |g| < 1 we will show the existence of Gq independently by deriving explicit expres- 
sions for it. In addition to integrable hierarchies these functions also appear in analytic 
number theory, see [60, Remark 3]. 

^-Multifactorials. In jl5] Nishizawa derives an explicit combinatorial formula for 
G q d ^ (see (1811) ). A more illuminating formula for our purposes can be expressed in 
terms of g-shifted multifactorials, q- multifactorials for short (Theorem H]). We streamline 
Nishizawa's approach by introducing them and systematically using generating functions. 

The classical g-shifted factorial as defined by Euler is (a; q)^ := n^o^ ~~ ac ?)- ^his 
notation is widely used in the theories of basic hyper geometric series [18J, modular forms 
and partitions [38]. A natural generalization essentially due to Appell is 

Definition 6 (^-multifactorials). A q-shifted d-factorial is 
(°j q)^ '■= 1 ~~ a ( no Q dependence); 

oo 

(a;q)<$:= ]J (1 - aq^"^) , \q\ < 1, d = 1, 2, . . . . (66) 

i!,...,i d =0 



For finite N we set 



( a . a )W.= (Q7) 



Our indexing convention is in line with [60J but differs from [46], our (a;g)oo is 
Nishizawa's (a;g)^ _1 ' ) . Also one should not confuse our ^-multifactorials with multiple 
q-factorials that have a separate variable q^ for each index. Our definition is recovered if 
all qk are set equal to q. We already used (a; q)oo to write partition functions in a closed 
form. 

Unlike the case d — 1 the finite version (a; q)^ is no longer a finite product of (1 — aq l ) 
unless d divides N. Nevertheless we have 

aw d=0 (i - ^ ii+ - +id ) n: =0 nx.., 4d _ 1= o(i - ^ 1+ '" +id ) 



_ ll^o( Q g';g)^ 1 _ TJ( nn i. n )(d-l) 



and by definition 



N-l 



M { ^=]l(a^ ] qYt 1 \ d>l. (68) 



i=0 
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A similar calculation gives another useful identity 

(«7;g)g = -^TT. (00) 



(a;q) { c d) 



Also (/-multifactorials can be expressed as a single product if we group the factors by 
powers of q. To this end it is convenient to introduce the Stirling polynomials 

that reduce to the binomial coefficients when z = N > n is a positive integer. Some of 
their properties are reviewed in the Appendix. From a generating function for Stirling 
polynomials ( IHHi) 



oo 00/ \ / 00 

E <» + - + " = E E 1 «"= E*' 

il,...,i<j=0 n=0 VuH \-id= n / \i=0 



n=0 v 7 



r. (71) 



In other words, x ) is the number of arrangements of d non-negative integers adding 
up to n. As a consequence we have single product representations 

00 00 

n (i-^ i+ - +jd )=n( i - a ^) rrl) ( 72 ) 

ii ,...,1^=0 n.=0 

00 

(ag5ff)SS = n(l-a9 ,i 

71=1 

In particular, l/(g;g)oo } = Y\n=A l ~ ?T® = m (q) 

as we claimed before. They also 
imply that — ln(a; q)^ is nothing other than the quantum polylogarithm of Kirillov 



d+n-2\ 
n-1 ) 



Ex. 2.5.8], [16] defined for \a\, \q\ < 1 as 

00 

--=Yl Kl \ k)s -v (73) 

Indeeed, by (JggJ 

- .n( ; ,)« = -E(' i + ^ 1 ) ln(l- = £ (" + ^ " M £ ^ 

n=0 ^ 7 n=0 ^ 7 fe=l 

OO OO ✓ , X OO 

= e t e ( r ) <«*>* = e t (i - ^ = Li « (ai «> <74 » 

fc=l n=0 ^ 7 fe=l 
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The identity holds for any a G C by analytic continuation. 

Closed formulas for G q . Typically when a classical object is g-deformed its theory 
becomes more complicated. Refreshingly, g-multigammas for \q\ < 1 are an exception: 
their theory is much simpler than its classical counterpart. The underlying reason is that 
it is possible to write finite products as ratios of infinite ones ( 167|) . the latter having a 
straightforward extension from integers iV to complex values z. An analogous attempt 
to write iV! = (jv+i)-(jv+2)-'(Ar+3)--- l ea ds to a nonsensical product of all natural numbers. 
The closest classical imitation is the Gauss product formula 



z 



T(z + 1) = lim -± — n 

that requires the rather involved theory of Weierstrass products. One can fruitfully turn 
things around and derive product formulas for T and classical higher multigammas 
from the g-deformed ones via a limiting procedure |47j . 

We leave as an exercise to the reader to iterate the functional equation for Gq and 
derive by induction from fl68l) that for non-negative integers N 

G«iN + 1) = fe 9 )S - ( ") („,)£> )£) 

=ri(^>» ,< ~ 1) ' + ' (/,) r'^is 1 -"'. (75) 

i=0 

Now N can be painlessly replaced with any complex z as long as we stipulate that 
q G D\R_. This way first the infinite products converge since q is inside the unit disk D, 
and second 

q z ._ e z\nq ig 

defined by choosing the principal branch of the logarithm. Of 
course, a priori there is no guarantee that the function so extended coincides with G q d ^ 
of Definition [51 We now prove that this is the case. 

Theorem 4 (Alternating formula). Let q G ©\R_ with D the open unit disk in C. Then 
Gq (z) is an entire function of z given by 

of (z + i) = n (?;?)£ ( ~ 1)i+1(dli) ( ' 1)d ; (76) 

i=0 

where (g;g)S are the q-multifactorials (1661) and f*) are i/ie Stirling polynomials ( 1701) . In 
particular, for quantum gamma and Barnes functions (d = 1,2) 

^ + teff)S? (7?) 

G,(* + i) = -^nr (g;g)« 2 (^tfg. (78) 



.RecaW t/iat (g; g)£? := 1 - g, (g; g)^ := (g; g)oo- 



35 



Proof. We have to verify that the righthand side of (I76p satisfies all three conditions of 
flo^]) . For the duration of the proof we use G { q d) only as an alias for this righthand side. 
Normalizations (i) follow by direct substitution of values. For the functional equation 
(ii) we compute 



i=0 i=0 



(^;« ( - 1)d (79) 



<f-l 



II (?; q) 



» (-D l+1 [( d i7-i)+( z ,:-)] (g;g) (d) (-1)^+1 [ 



OO 



(-!)« 



i=0 



(q z ;q) 



(d-i) 



By ( 1901) with n = d — i the sum in brackets is just (/■) so the first two factors combine 

it f 

(q z ; q) l £ 



into Yli=o (?i 9)<» • For the last factor, (16"9"j) with a = q z yields 



(q z ;q) 



(d-i) 

00 



so fl79|) becomes 

G^\ z )Gf{z) = f[(q;q) { 2 { ~ iy+1{dU \q z+1 ;q) { £ = Gf(z + 1). 

It remains to verify the log-posit ivity condition (iii). Taking the logarithm we have 



in Gf>{z + 1) = (-!) m Hq; ?)2 (/ •) + (-!)" Hq z+1 ; q)£ ] - 



(80) 



Since (/) is a polynomial in z of degree d — % the sum above is a polynomial of degree 
d and its <i + 1-st derivative vanishes. Recalling (IT^j) one obtains 

' \nGf\z) = (-1)^—— Li d+1 (<f;g). 



dz d + l 

d df 
It is helpful to notice that -rf(q z ) = \nq x — 

dz dx 



dz d+1 



and iterating 



x=q* 



Ad+1 



lnq) d+1 ( x-f - ) f{x 



d 



d+l 



Now for \x\ < 1 by (J73 

d+l 



{ X Tx) LWi(x;g) = g n(1 ^ n)d (x 



f/.r 



d+l 



x=q^ 



Ell 
nil - a' 1 



n=l 



q n ) c 



x 
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and therefore 



^TT ^GfO(z) = (-lng)^ ( — ) <T for |g« 

n=l 



< 1. 



If < q < 1 and z is real positive then — In q > and < g z < 1 so the last expression 
is positive by inspection. □ 

The structure of G q is perhaps most transparent in flHUl) . We have the main term 
(— l) d ln(g z+1 ; that depends on a = q z only and the anomaly term polynomial in z 
of degree al. For d = 2 we recognize ln((f + 1 ; g)£ } as exactly the reduced free energy of 
the resolved conifold, see (1531) . The main/anomaly structure of Gq becomes even more 
transparent if we introduce a generating function for the entire hierarchy 

oo oo oo 

in Gf{z + i)t d = (i+ ty LWi(g; q) t d - (-!) d LWi(g^ +1 ; ?) ^ 

One can check that the sum converges for \q\, \q z \ < 1, |t| < 1 — This formula displays 
very clearly the alternating nature of G q and the distinction between the main and the 
anomaly terms. The dependence of the anomaly on z is regulated by a universal term 
(1 + t) z that explains the convolution structure of the anomaly sum in (180]) . 

This structure is quite common for functions appearing in quantum field theory. The 
main terms usually reflect the expected symmetry of a system while anomalies require 
additional choices to be defined. For instance, Witten's Chern-Simons path integral is 
not a 3-manifold invariant because of the framing anomaly [H EI] but can be turned 
into one by choosing the canonical 2- framing to cancel the anomaly [58] • In our case the 
obvious extra choice is that of a logarithm branch to define z := log a. Disregarding 
anomalies comes at a price: fixing the canonical 2-framing would complicate the gluing 
rules and Gq would not obey a simple functional equation. The Reshetikhin-Turaev 
normalization may be to blame for nasty prefactors in flMl) that spoil the duality. It 
might be of interest to find an analog of 2-framings on the Gromov-Witten side and 
compare answers without making choices, even canonical ones. 

From the alternating formula (!76j) we now derive a single product formula for Gq® 
originally given by Nishizawa |45j. 

Theorem 5 (Nishizawa product). Let q G D\M_ then quantum multigammas are given 
by 



z — n\ 



GfHilHl-f !]- ^T^T ford>l. (81 



n=l {l-q***)\4-V 
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In particular, for quantum gamma and Barnes functions (d = 1,2) 

r 9 (z + i) = (i-g)-n t^t (82) 

n=i ^ y / 

G„(* + 1) = (1 - g)"^ J] n_lJ + n (83) 



n=l 

Proof. Applying (!72|) to the first product in (176|) we get 



..^«> <-i>'+'U.) 



00 

= ]J(l-q n+1 )^- 1)t+i rr )(/-,). (84) 



i=0 i=0 n=0 

00 



71=0 



As shown in the Appendix (19T1) . when d > 1 



i=0 



i^ + i-1^ f * \ = f(^V)> n>l 



Splitting the n = and n > 1 factors in (1841) we see that it equals 

00 00 
(l- g )(3=i)"(S) n(l-g" +1 )( Z -i 1 ) = (l-g)-© n(l-g" +1 

n=l n=0 

= (l- 9 )-Cy JJ(1 - g»)C3 

71=1 

The second factor in (1761) can be transformed into a single product using (172)) again 

00 00 

_ ^(-^(Tr 2 ) = JJ(i-^ 



z — ra \ 

-1) 



n=l n=l 



where we also applied (J87j) and ([50]) to get the second expression. Multiplying the right- 
hand sides of the last two formulas gives the claim. □ 

Formula (l82l) is the standard expression given for the quantum gamma function, see 
e.g. [5], and (153]) is its closest quantum Barnes analog. Although Nishizawa product (18"TT) 
may appear prettier than the alternating formula (I75|) it completely hides the graded 
structure of g-multigammas. 

What about |g| = 1? By Theorem H]g-multigammas are entire functions of z for \q\ < 1. 
Things change if we venture onto the unit circle. As one can judge by the example of 
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the classical Euler's T where q — 1, (jHSD still defines a unique function of z but this time 
it is only meromorphic. The appearence of poles prevents infinite products (181 p from 
converging. When |g| = 1 but q is not a root of unity G q = T q is constructed explicitly 
in [1H] via Shintani's double sine function (see [60]). This case is complementary to the 
classical one q = 1 and the poles are located at the points — n — m/r with n, m e Z and 
q = e 2mT . Clearly, it is desirable to have a similar construction for = G q especially 
when q is a root of unity. Indeed, this case is most relevant to the Chern-Simons theory 
where G q is essentially the partition function of 5* 3 by Theorem [3] and probably appears 
as a factor in partition functions of other manifolds. 

There is a small window through which we can peak at what happens on the unit 
circle for any q. When z = N is a non-negative integer the product (IHTj) terminates. 
In the context of Chern-Simons theory N is the rank of sljvC and we already saw this 
phenomenon for the quantum Barnes function in Lemma [TJ This also generalizes the fact 
that although T(z + 1) can not be simply expressed as an infinite product for complex z 
we have nonetheless T(N +1) = 1 • • • N. 

Corollary 2 (Nishizawa). For a non-negative integer N and any q G C one has 

N 

G q d \N + !) = (!- Yl (1 - g n )^ : "). (85) 



71=1 



Proof. By the Nishizawa product for Gq 
{1-qp) G q d \N + l) 



rr=i (i - i n 



N-n\ 
d-1 ) 



rELi(i-? 



n+N 



N-n\ jy 

d_1 ' T— T . _ . (N-n\ 

d-1 ) 



IK 1 -"' 



71=1 



This gives a proof for \q\ < 1 but since the righthand side is a polynomial in q the general 
case follows by analytic continuation. □ 

Note that by (JS5J G q d) (N + 1) = for N > n if q is an n-th root of unity but this 
never affects Chern-Simons invariants since for them n = k + N with positive integer k. 
This termination phenomenon sheds some light on why so far Chern-Simons invariants 
have only been defined for integral z = N, when it can be interpreted as the rank of a 
quantum group or an affine Lie algebra. To include complex z one should probably study 
more involved algebraic/analytic structures. Reversing the perspective, we observe that 
termination is what makes G q meaningful in Chern-Simons theory that in its current 

form produces values only at roots of unity. Note that with the exception of (1 — g)~( d ) 

all anomaly terms in ( ITS]) are required for termination. In particular, the ratio ^ aq,q, ^ ) 

with q z = a that we recover from Gromov-Witten/Donaldson-Thomas theory, does not 
suffice for duality to even make sense. It does not terminate for z = N and becomes 
meaningless for q = e 2m /( k + N ) ) which is where the Chern-Simons invariants are defined. 
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Conclusions 



Recently the phenomenon of holography has become prominent in physics [57] . Roughly, 
the idea is that to every theory on a bulk space there corresponds an equivalent theory 
living on its boundary. The most famous example is the AdS/CFT correspondence of 
Maldacena but one can certainly trace the analogy back to the classical potential theory, 
where a harmonic function is recovered from its boundary values. We see a toy example 
of holography playing out on the ranges of Calabi-Yau invariants. The master-invariant is 
a holomorphic function in the bulk (the unit disk) and Donaldson-Thomas theory comes 
closest to being a bulk theory by giving its Taylor coefficients at 0. Chern-Simons theory 
gives its values on the boundary (the unit circle) and Gromov-Witten theory also lives 
on the boundary but via asymptotic coefficients at 1. Recall that q = e tx and in terms of 
the string coupling constant x we are dealing with periodic holomorphic functions on the 
upper half-plane. Their Fourier coefficients are given by Donaldson-Thomas invariants 
and values at the cusps of SL<{L are given by Chern-Simons invariants. This brings to 
mind the classical modular forms [B] and indeed the relationship between them and the 
corresponding boundary objects was interpreted recently as an example of the holographic 
correspondence [39J. Although our functions are not modular in the classical sense they 
do have some modular transformation properties |34j . This links large N duality to a 
well-known problem in analytic number theory and perhaps this can aid in its proof. 

There are some serious difficulties to be resolved. First, values of a holomorphic 
function at roots of unity are not enough to recover it uniquely. There are plenty of 
holomorphic functions on the unit disk that vanish at all roots of unity. Take a modu- 
lar cusp-form / on the upper half-plane [6] for example and consider f(lnq/2iri). The 
same is true of an asymptotic expansion at a point on the unit circle or even a collec- 
tion of asymptotic expansions at every root of unity [54] . In other words, absent extra 
data Chern-Simons and Gromov-Witten theories do not recover the master- invariant. 
Donaldson-Thomas theory has a converse problem: it sure gives a function on the unit 
disk but one that is singular at each point of the unit circle. For the resolved conifold 
we reconciled the results by 'completing a pattern' by hand but this will not work in 
general. A large part of the problem in proving equality of partition functions is that as 
mathematically defined they are not quite equal. 

We contrived to make large N duality work for the resolved conifold. What about 
other cases? Similarities between the topological vertex and the Reshetikhin-Turaev 
calculus outlined in Sections I3H41 are promising but there is a vast difference between 
toric webs and link diagrams representing the dual objects. This ought to be expected 
of course since not all large N duals to 3-manifolds are toric and certainly most toric 
Calabi-Yaus are not dual to any 3- manifold. However, not all is lost. Link invariants 
are known to extend to invariants of knotted trivalent graphs [44J and these include toric 
graphs that are also trivalent but planar. Then there is the case of lens spaces (cyclic 
quotients of S 3 [52J) where the duals are known and toric [21 [23]. Moreover, there exist 
physical generalizations of the topological vertex to non-toric threefolds [32] some of 
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which are dual to non-cyclic spherical quotients. 

Even if diagrammatic presentations are reconciled there is yet another hurdle to 
wrestle with. Topological vertex expressions are infinite sums over all partitions V. 
Reshetikhin-Turaev expressions are finite sums over partitions in the (N — 1) x k rect- 
angle "P^r-i an d a d hoc tricks are required to turn them into functions of q = e 2m /( k + N ). 
Note that the Reshetikhin-Turaev building blocks W Xfl (N;q) are already expressed in a 
desirable form. It is when the gluing rules are applied that q has to be specialized to 
a root of unity to render the sums finite. Analogous infinite sums naively diverge but 
so would sums of Wx^q) in the topological vertex without the convergence factors like 
(— a)' A ' in (I52I) . After the sum is performed it is a again possible to turn it into a function 
of N, q up to framing and volume factors. Recall that values at roots of unity do not in 
themselves determine a holomorphic function on the unit disk appearing in the dualities, 
and Donaldson-Thomas theory as a bulk theory is incomplete. An enticing possibility is 
that Chern-Simons theory itself can be turned into a bulk theory. 

Conjecture: There exists a universal Chern-Simons TQFT that assigns holo- 
morphic functions of q, z to links and 3-manifolds with q G D. These functions 
extend to roots of unity in q and specialize to ordinary Chern-Simons invari- 
ants at rank iV and level k upon specializing to z = N, q = e 2m /( k + N ) U p to 
normalization factors. Partitions serve as colors of link components and the 
gluing rules are expressed via sums over all partitions. Donaldson-Thomas 
invariants of the dual Calabi-Yau threefolds are Taylor coefficients at q = 
of these functions and Gromov-Witten invariants are their asymptotic coeffi- 
cients at q = 1. 

This conjecture formalizes the idea that extra extension data for invariants comes from 
algebraic restrictions required by the axioms of TQFT [7J. A corollary of this conjecture 
that seems to be within reach is that sums of the type 

w XXl w XlM ---w Xn , 

Ai,...,An 

can be represented as boundary values of sums 

w XXl w XlX2 ---w Xn , r i Ai i + - + iH 

\i,...,x n ev 

where r is a convergence factor. These and similar sums appear in Chern-Simons in- 
variants of Seifert-fibered 3-manifolds [4*0l [53] and such representation would facilitate a 
proof of large N duality for them. 
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Appendix: Stirling polynomials 

The Stirling polynomials [38], 1.2.11 are defined as 

"V\ / z\ z(z — 1) • • • {z — n + 1) 



0/ \n j n\ 



n > 1, 



and reduce to the binomial coefficients when z = N > n is a positive integer. Their 
generating function is 



n=0 v 7 



56) 



In agreement with (186]) we always assume (^) = for n < or n > N. Negative 2 is 
also allowed 

- A '.ir( z+n - v ). (87) 



n J \ n 

Combining with the generating function (1861) we have 

CO ✓ 1 

z + n — 1 



n=0 



The coefficients s(n, k) in 

n.. k.} . f 7. 4- r> — 1 \ . .sfn. 

z fe (89) 



fz\ = j2 s ( n > k ) z k. fz + n-l\ = j^^n-k 4n, k) j, 
fc=o n ' \ n J k=0 7 



are known as the Stirling numbers of the first kind [55] . 1.2.11 . The following identity is 
easily derived via generating functions 



K n — 1 J \nj \ n 
For the convenience of the reader, we prove a harder one needed in Theorem El namely 

6f V « )\d-i) n = o. (J1) 

Proof. Consider the generating function 



fv(£(-i)- + tr%^H t,+ ' 

d=0 \i=0 v 7 v 7 / i,j=0 



n 



E( _ 1)<+ . " + -i ,. E Qw, 

i=0 v 7 j=0 v/ 



oo 
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The second factor is obviously (1 + t) z by (15oT) . Assuming n > so that ( n ) = the 
first factor can be rewritten as 



( £ /,-! + („ +1,-1 )( _ t) ,- 1=t(1+t) 



-n-l 



i - 1 

i=i v 

by (1551) . Multiplying them and expanding in the powers of t 



We may conclude for n > 1: 



n + z-A / z \ = f(— d>l 
n J \d-i) |0, d = 0. 



E(-W 

2=0 ^ 

The case n = is special since = 1 by definition. Therefore, 

i=0 ^ ' 

and expanding as above we obtain for n = and all d 



By assumption we have d > 1 so applying (190|) 



□ 
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